Coherent control of artificial molecules using an Aharonov-Bohm magnetic flux 
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Bonding and anti-bonding states of artificial molecules have been realized in experiments by di- 
rectly coupling two quantum dots. Without a direct coupling between two nearby quantum dots, 
here we show that a continuous crossover, from symmetric to anti-symmetric molecular state, can 
be achieved by changing the flux through a double quantum dot Aharonov-Bohm (AB) interferom- 
eter. We explicitly present the flux- dependent real-time processes of molecular-state formation. In 
contrast to the transport current, which has a 2tt period, the quantum state of the DQD molecule 
has a 4:71 period in the AB flux. 

PACS numbers: 73.63.Kv, 03.65.Wj 



It is important to tailor quantum states, especially, 
to control the coherent phase between two superposition 
states. In the past decades, artificial atoms and molecules 
in solid-state systems, such as superconducting Joseph- 
son junctions pQ and semiconductor quantum dots (QDs) 
[2]-[4] have provided novel platforms for exploring such 
quantum-coherent effects. Due to the tunability of vari- 
ous electronic couplings, double quantum dot (DQD) sys- 
tems, which are archetypes of artificial molecules, have 
attracted considerable attention. Using direct-tunnel 
coupling, the coherence of charge states have been in- 
vestigated with Aharonov-Bohm (AB) interferometers in 
recent experiments [5j[6]. However, the couplings to the 
electron reservoirs (electrodes used for measurements and 
controls), generally induces decoherence to the quantum 
state of the DQD molecule. Here, we show that for an 
uncoupled DQD in an AB interferometer, such decoher- 
ence can be suppressed with an asymmetric design of the 
device geometry. Therefore, by solely tuning the AB flux, 
the coherent control of the DQD molecule (from the sym- 
metric to the anti-symmetric state) can be realized. Fur- 
thermore, we find that the period of the quantum state 
of the DQD in the AB flux is 4tt. The transport current, 
obtained by averaging the DQD states, possesses a pe- 
riod of 2tt. The coherence of the DQD molecular state 
and the coherence of electron transport therefore mani- 
fest themselves fundamentally different through the AB 
flux. 

Coupled DQDs have been theoretically proposed as 
qubits pj E] and experimentally realized [3] 2J I9HT2]. 
The feasibility of realizing various one- and two-electron 
molecular states with tunable tunneling and exchange 
couplings was demonstrated [T3j [14] . Furthermore, the 
coherence of electron transport has been investigated 
with a single QD in AB interferometers [15]. Combin- 
ing an inter-dot tunnel coupling with a magnetic flux 
has also been studied theoretically [T6UT8] and exper- 
imentally [5j [6] [19] [20]. In particular, controlling the 
molecular-state through AB phases is of recent exper- 



imental interest [5 . Although tunneling to the elec- 
trodes may be turned off to avoid the electron-reservoirs- 
induced decoherence, such tunnelings are indispensable 
for the AB effect. Thus controlling the molecular coher- 
ence through the AB flux is a new challenge. 

Here we consider uncoupled DQDs embedded in an 
AB interferometer, as sketched in Fig. [T] In contrast 
to previous theoretical studies, which focus on quantum 
transport [17] [18] , here we directly exploit the quantum 
state of the artificial molecule. By explicitly analyzing 
the decoherence through the AB flux, we deduce the 
proper geometry of the DQD for coherent control over 
the molecular states. The time-resolved formation pro- 
cesses of various molecular states, determined by different 
AB fluxes, are explicitly presented. 




FIG. 1: (color online). A schematic diagram of a pair of un- 
coupled quantum dots in an Aharonov-Bohm interferometer. 



The model system and its exact solution. — To focus on 
the influence of the AB flux on the quantum state of 
the artificial molecule, we consider only polarized non- 
interacting electrons. The total Hamiltonian of the sys- 
tem is conventionally [21 given by H = H s + He + Ht, 
in which 7i s = Eia\a i describes an uncoupled DQD 

and 7~L E = ^Z a k e ak c ak c ak is tne Hamiltonian for the 
leads with a = L(R) labeling the source (drain) lead, 
and 1-Lt = Yjj a k[^j^ c ak a j + H.c] depicts the coupling 
between the central dot system and the leads. Here a\ 
(aj and c* ak (c a k) are the electron creation (annihila- 
tion) operators for the electronic levels i and k in the 
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dot system and the lead a, respectively. The tunnel- 
ing amplitudes harbor the applied magnetic flux <fr via 



V* — 
V 1L — 



V 2L = \V L \e^/\ and V 1R = 



V 2R — 



\V R \e 



i<f>/4 



where (j) = 2tt^/^q and $0 = h/e is the flux quantum. 
The line-widths induced by tunneling are then given by 
T a = 27r\Vj a \ 2 Qa, where g a is the density of states in the 
lead a. The DQD molecular states are governed by the 
following master equation [22] : 



dt 



p(t) = -i[U s ,p(t)} + + £Z>V)W)> (!) 



with / being an identity matrix and 



(4) 



where £f a (t) are the superoperators describing dissipa- 
tions and fluctuations induced by the tunnel coupling to 
the electrodes (for details, see Ref. [22]). Denoting the 
state of the empty DQD by |0), one electron on the first 
and the second dot by |1) and |2), respectively, and the 
state of both dots occupied by |3), the density matrix 
p(t) can be generally expressed as 



pit) 



fPoo(t) 



V 



\ 

Pu(t) pi 2 (t) 

P2l(t) M*) 

p 33 (t)J 



(2) 



where pij = (i\p\j) with z,j = 0,1,2,3. The molecu- 
lar state, featured as one electron in the DQD shared 
between the two orbit als of the dots, is embedded in 
the central 2x2 block matrix of Eq. ([2|. In particu- 
lar, the coherence between the two atomic orbitals of the 
DQD molecule is characterized by the off-diagonal ele- 
ment, P21. To see how molecular states in this DQD are 
formed in time, we solve the master equation with the 
initial preparation of empty DQD, namely, poo(0) = 1 
and pij(0) = 0, for all i 7^ 0, j 7^ 0. The explicit solution 
of each matrix element gives 

pu(t) = vu(t) - detu(t), 
P22O) = v 22 (t) - detu(t), 

Pl2(t) = Vi2(t), p2l(t) = V 2 l(t), 

p 00 (t) = det[I - v(t)], p 3 s(t) = detv(t), (3) 



is a 2 x 2 hermitian matrix, where f a (w) is the Fermi 
distribution function of the reservoirs, the upper (lower) 
sign is for a = L (R), and u(t,u) = J tQ dTe lUJ ^~ T ^u(T) 
with 



u{t) = exp 



r c *fa) 



iE 2 



(5) 



Here we have defined T c (<p) =T cos(0/2) + iST sm((j)/2) 
with r=r^+r^ and ST=Tl— Tr. The functions u(t) and 
v(t) are indeed the retarded and correlation Green func- 
tions in the Schwinger-Keldysh nonequilibrium Green 
function theory [23]. The AB flux 0, the coupling asym- 
metry (5r, the non-degeneracy SE = E\ — and the 
nonequilibrium dynamics from the electron tunnelings, 
all influence the consequent quantum states of the DQD 
molecule. 



Coherent phases controlled by the AB flux. — To have 
a clear picture of the coherence of the DQD molecular 
state, let us look at the off-diagonal matrix element pi2(t) 
in Eq. ^ first, in the steady-state limit (t > T" 1 ). The 
general solution gives 
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where j ((/)) = yj T 2 cos 2 (0/2) + ST 2 sin 2 (0/2) - SE 2 and 
r±(0) = 2 _1 (r ± 7(0)). Here, we have also applied a 
bias Pl = eV/2 = —p R at zero temperature. The full 



complexity of decoherence is revealed through Eq. (|6]). 
Due to the severe decoherence in such system, it can be 
proven [24] that the coherent phase ip (in the off-diagonal 
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matrix element p2i = \p2i\e Z(p ) between the two atomic 
orbitals can only take the values of 0, ±7r/2 or ir for ar- 
bitrary flux. This applies for the often-used condition of 
degeneracy SE = and symmetric coupling ST = 0. The 
decoherence-induced discretization of the coherent phase 
hinders the manipulation of the coherent phase of molec- 
ular states. However, when the DQD is non-degenerate 
and couples asymmetrically to the left and the right leads 
(ST 7^ 0), we find that the coherent phase <p can be con- 
tinuously tuned by the AB flux. 

In order to achieve typical molecular states, the DQD is 
set at degeneracy (SE = 0). The second term in Eq. (|6| 
vanishes. Equations ( |3|6[ ) show that the formation of 
molecular states is essentially determined by the applied 
bias and the coupling asymmetry to the source and the 
drain. The basic setup of zero bias (which is used for 
examining quantum transport) leads to p2i =0 (because 
the DQD is in equilibration with the reservoirs) and is 
not interested here. With a large bias, eV ^> T, we find 

f | [(*r/r) cos * -isinf] if <t> + 0, 
P21 = <^ " " (7) 

[ !(i + *r/r) if = 0. 

Equation ([7]) clearly shows the controllability of the co- 
herent phase between the two atomic orbitals of the DQD 
molecule through the AB flux. It also explicitly reveals 
the necessity of the asymmetry in couplings, ST ^ 0. In 
the case of symmetric coupling, ST = 0, Eq. ^ shows 
that the real part vanishes for ^ so that the coher- 
ent phase if is localized at 7r/2, except for <\> = where 
the coherence phase is restricted to 0, as it has been 
pointed out in [24]. With the larger asymmetry ST, the 
coherence amplitude |p2i| linearly increases and the co- 
herence phase is continually driven by the AB flux, as 
seen from Eq. ([7| [25]. By setting ST < T, we obtain 
P ~ where 

1^(0)) = ^=[|1) +exp(-i0/2)|2)]. (8) 

A continuous transition from the symmetric, |^(0)) = 
(|1) + |2)) /a/2, to the anti-symmetric state, | < 0(d=27r)) = 
(|1) — |2)) /a/2, is achieved by changing the AB flux, as 
shown in Fig. [2] (see captions) . Interestingly we find that 
the period of the state of the DQD molecule is 47r, rather 
than 2tt in the AB flux as one expected. 

Real-time processes of molecular- state formations. — 
The full information of the quantum state of the DQD 
molecule at finite temperature is depicted by the time- 
dependent reduced density matrix. We can write the 
central block matrix of Eq. Q as 



p q (t) = -[l + r(t)-a] 



2 M*)+P33 (*)]/, (9) 



where a — (<r x , a y , a z ) consists of the Pauli matrices and 
r(t) = 2{Rep2i(t)Jmp 2 i(t),p 11 (t) - P22O)}, is the po- 
larization vector for the molecular states. So the dy- 
namics of molecular-state formations can be visualized 
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FIG. 2: (color online). Control of the coherent phase of the 
DQD molecule by the AB flux. The explicit solution of p2i in 
the steady-state limit is shown by the "blue diamonds" on the 
central panel. Each diamond corresponds to an AB flux value, 
taken from <j> — to <j> — ±2tt with tt/8 steps. The wavefunc- 
tions on the DQD molecules are illustrated for various values 
of the AB flux. A DQD is indicated by two circles with cen- 
ters connected by a dashed line (no inter-dot coupling) below 
the diagrams for cj) = —3tv/2. Both the real (red) and the 
imaginary (blue) parts are shown, so one sees how the AB 
flux changes the coherent phase between the atomic orbitals. 
Other parameters are SE = 0, eV = 6T at ksT = T/20, which 
are also used in the following figures, unless specified. 



through the motion of the polarization vector with the 
Bloch sphere. Also, the leakage out of the one-electron 
state-space can be easily seen from the term proportional 
to the probability of the empty and the double occupied 
states, p 00 (t) +P33(*)- 

In Fig. [3] we plot the evolution of the full reduced 
density matrix of the DQD molecule. Initially, the DQD 
is prepared in an empty state, poo(0) = 1 as shown by 
Fig. [3jbl,cl,dl) and r(0) = given in Fig. [3jal) (where 
the length of the red strip is zero). After injecting elec- 
trons from the left and the right reservoirs, poo decreases 
[see Fig.|3jbl) to (b3)] while the electron occupation and 
coherence increase with time [see plots (bl) to (b3), (cl) 
to (c3), and also (dl) to (d3) in Fig. |3j. The coher- 
ent phase (f between the atomic orbitals has been fixed 
shortly after the electron injection [see Fig. j3^al) to (a4) 
and also (d2) to (d4)]. Then \r(t)\ grows in time with 
fixed (/?, and finally a stable molecular state, p w l^)^, 
where = (|1) + e -^/ 2 |2)) /a/2, is reached in a short 
time, of about 3r _1 . Note that due to possible leak- 
age, see Fig. [3^b3) where poo has a small finite value, the 
DQD is not in a perfect pure state. But the situation 
can be optimized by changing the bias and the coupling 
asymmetry, as shown by Eq. ([7]). 

To better understand the role played by the AB flux, 
we show the time evolutions of p2i in Fig.[4]under various 
values of (f). Figure [4jal,a2) show the process of coher- 
ence generation for a strong asymmetric coupling (with 
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r=0.6/r 




FIG. 3: (color online). Typical process for forming molecular states. The dashed black line in (al) to (a4) is the trajectory 
taken by r(t) from t = in (al), starting from the origin, to t — 3/T in (a4), where it almost touches the surface of the 
sphere. The red strip in each plot is the trajectory up to the corresponding time points, as shown above the spheres. From 
the trajectory, we see the coherent phase ip (which is the angle made by r(t) with r x axis) has been fixed after the electron 
is injected into the DQD. Plots (bl) to (b4) display the real part of the reduced density matrix of the DQD system, while 
the imaginary part is plotted in (cl) through (c4). The coherent phase ip between the two atomic orbitals is better visualized 
through the vector plots (dl) to (d4). Every arrow represents an element of the reduced density matrix pij, whose horizontal 
projection stands for the real part and the vertical projection stands for the imaginary part. The AB flux here is = — 7r/2. 



ST = 0.9T). The rate of approaching steady-coherent- 
molecular states is only weakly dependent on the flux. 
The stable molecular states are soon reached after a few 
This is totally different from the symmetric cou- 
pling (ST = 0), see Fig. gbl,b2). It shows the flux- 
dependent decays of Rep2i, due to the severe decoherence 
in the symmetric coupling. Therefore, the coupling asym- 
metry can strongly suppress the decoherence induced by 
electron tunnelings, and make the coherence control of 
the QDQ molecule feasible. 



Discussions. — The general solution shows that the 
quantum state of the DQD molecule has a period of 
47r in the AB flux. It is an intrinsic property of this 
pseudo-spin system, independent of the coupling geom- 
etry and the bias configurations. Besides, we have also 
calculated the tunneling current to reservoir a = L,R 
within the same framework [22] with the result I a (t) = 
e tv s [Cf a (t)p(t)]. The steady-state transport current 
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FIG. 4: (color online). The time evolutions of p2i. The plots 
(al,a2) for 5T = 0.9T and (bl,b2) for 5T = 0. Also (al,bl) 
give Rep2i, and (a2,b2) give Imp 2 i- 
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I = \{I L - I R ) is then given by 

where the transmission coefficient is given by 



T(w,< 



(r 2 -,5r 2 )[ w 2 cos 2 



iffisin 2 f] 



(10) 



(11) 



By taking ST = it reproduces the result in Ref. [21] . 



Equation (11) clearly shows that the transport current 
has a period in the AB flux of 2tt. This 2tt period, as 
a feature for the coherence of transport, is well known 
and has been observed in experiments [15j [19j [20] . The 
47r period, a nontrivial character of the quantum state of 
the DQD molecule, requires further experimental inves- 
tigation. Note that although the coherent phase of the 
off-diagonal density matrix element is gauge-dependent, 
the AB flux dependence of the coherence phase and its 
periodicity are both independent of the gauge choice. 

In summary, we have demonstrated the effectiveness 
of the AB flux for the coherence control of DQD artifi- 
cial molecules. We have analyzed the AB flux-dependent 
conference controlling, through the asymmetric coupling 
of the DQD to the electron reservoirs. When a large bias 
is applied with a strong asymmetry in couplings to the 
source and the drain, coherent control by the AB flux can 
be easily achieved. The decoherence induced by the elec- 
tron tunnelings can be efficiently suppressed. We also 
find that the period of the quantum state of the DQD 
molecule in the AB flux is 47r. The revelation of the 
underlying quantum-coherence of the molecular states is 
thus beyond the usual transport measurement. The veri- 
fications of these molecular states would rely on a suitable 
quantum-state-tomography protocol for further investi- 
gations. We hope that this theory for artificial molecules 
could inspire new experiments on coherence control of 
molecular states via the AB flux, and become also use- 
ful for the quantum emulation [26 of artificial molecular 
processes. 
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